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Piiring the sixth quarter of this project, the DST liidicer have been com- 
piled for the time period that Magsat was in orbit, and magnet ic tapes con- 
taining these indices have been delivered to both tlie tiala c»*nter and the 
Magsat project manager. 

As an alternative method for verifying the DST correctional scheme pro- 
posed earlier In tliis Investigation, another appr»)acl> has been developed for 
separating the portions of the magnet^!’ field measured by Magsat tlwit arise 
from internal and external sources. 

In this approach we begin by considering a set of functions wliich satisfy 
Dirichlet's problem for the circle: 


u (r,6) 




cos n0 + b 

n 


sin nU) 


where the source is Inside the circle, and 
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for sources outside the circle. 

These functions are harmonic and tell us that If we know the value of a 
potential function on a circular boundary, we know Its value everywhere inside 
and outside the circle. The two equations describe the cases where the sources 
are internal and external. 

In the case of the Magsat data, measurements are made both above and below 
a circular boundary due to the elliptlcity of the satellites orbit. Therefore, 
impliciL in these measurements is Information on the attenuation of each fre- 
quency in the magnetic field for a particular orI>it as .1 function of altitude. 
Adcscription Is given Iiere of an attempt to s«-paiMti- the inti-rnal ami external 
fields measured bv an orbiting satellite. 

A particular measurement of the earths magnetic field at a distance of r 
kilometers from tlie center of I lie earth can be expressed as 



V = 


s I n 
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To test this method, a set of sample cuef f ie lent •. witc usnl to <-oiii|uiti' 
field values alotiK a simulated satellite orbit. This data was I lu>n used to (tv 
ti> ree»>ver the ori);inal eot*f f ie lisit s used to >^eiieia(c Clu’ I ieltl valiiets. The 
satellite orbit was given by r - 6800 + 100 sinO where 0 is the central angle of 
the orbit and max and min values of r were 6700 km and 6900 km with 6800 km the 
mean value. To compute the field 6 harmonics were used for the inner terms and 
six for the outer. No constant term was used for tlie input but one was computed 
for the output. The input coefficients were solved for in two ways: 

(1) liA-ersion of a 25 x 25 matrix 

(2) A modified version of the least squares solution wliich utillxed the matrix 
inversion lemma, to obtain the recursive least squares solutfmi. I'or this 
case, 2048 points were used for input. Note that in using this nu*thod it 
... iu)t necessary to invert a matrix. 

Table 1 shows the results of solving for input coefficients using r(lnt) = 
r(ext) = (»800 kni. It can l>e seen tliat tiie matrix inversion metiiod gives the <-or- 
rect values for tlie 6tlt liarmonic, hut tlie lower order terms heconie increasingly 
bad. the recursive least squares method doesn't do as well, Init is in the right 
ball park tin* tin* sixtii liarmonic. 

In tables 2 tlirougli 4, the same input data was used, but the radii of the 
circles on which tlie internal and external boundaries were defined was varied. 

In table 2, r(int) = 6400, (approximately 1 eartli radii) and r(ext) = 5 x 6400 
(approximately 5 earth radii). In the matrix inversion solution the external 
terms blow up, while only the 5th and 6th terms are approximately right for the 
Internal coefficients. The recursive least squares solution gets all the inter- 
nal terms correct (which are all that really contribute anything to the field 
at the satellite orbit in this case) but the external terms never really get off 
the ground. This seems to be a peculiar property of the recursive least squares 
method. It is probably due to the fact that the very small attenuation coef- 
ficients for the external terms tend to suppress their contribution completely. 

riu“ values in table 3 used 6400 ftir r(int) ami 7:’00 for rfext). 'I'luis t tie 
Internal .ind external generating lioundaries were approximately equidistant from 
the satellite orbit. Again, tlie matrix inversion was only rig.lit for tlie iiiglu-r 
order terms, wliile the recursive li-ast scpiares seemed to split uji tlie coefficient 
approx iniate 1 V i veiilv Ik tween internal and external tc'rms. 


If wo t.iko a series of 4N+1 measurements at equal angular Increments tlirouKli 
the courst' of one orbit • (1*1,2,. . .,4N-fl), we )»i*t 4N+I e(|uat ions and nn 

equal number of unknowns. These equations can In- written In matrix notation 

. 1 ;. w!,i r, , 


y ■ Mx 

y “ 5'i 
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X - (aQ, ai* bi 
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If there are more measurements than unknowns, a least squares solution can be 
obtained using the relation 

'' T -1 T 
X - (M M) M Y 

where y is an mxl column vector and m>4N+l. 


It should be noted that this method is similar to that used by Henderson & 

Cordell in adjust Inn in>tentlal field data to a eommoii altitude. The <liffereiuM' 

beinn that liere the data is considered for the prohli*m of t Ite circled rather 

than for a strainht line or plane; i.e., the attenuation factor is of tlic* form 
r n 

frCi) rallur than e . In anv case, this method lan be used to adjust satel- 

lite dat.i to a common altitude, and, if it is appl iid to ci>mpoiu‘nt d.ita r.itlier 
lian L : a 1 ;i< 'd viata, t lu assumption of harmoni«'itv is valid. 


laMi? A list'd r(int) ■ 6400 km and r(ext) • 2 x r(lnt) 

I'aMe 5 was the same with different Input eoef f le lenls. 

In table 6 r ( lnt)“ r (ext)" 6800, as in the first case. However, here the 
equation for the satellite orbit was changed to 

r - 6800 + 3200 x sin 0 

In this case, the last 4 harmonics came out correctly In tlie matrix Inversion, 
suggesting that the matrices are not very well conditioned for the case where the 
Magsat orbit was approximated. The recursive least squares also performed some- 
what better. 

Several other methods are available to solve for these coefficients and 
will be reported on shortly. 
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Internal 


External 


Harm 

Cos 

Sin 

Cos 

Sin 

1 

- 35.09 

22.04 

- 2.92 

“ 0.56/ Input 

2 

22.28 

4.44 

1.86 

- 5.48 1 Coef f 

3 

- 22.95 

26.89 

2.29 

11.16 ' 

4 

18.87 

0.20 

- 2.61 

- 7.75 

5 

7.86 

27.34 

1.74 

6.68 

6 

- 0.69 

17.99 

- 1.50 

- 0.61 


.) 


a - -0.25 
o 


1 

- 19.07 

1 ! .05 

- 19. OK 

10.78 

2 

11.83 

- (J.51 

12.13 

- 0.83 

i 

- 10.28 

19.06 

- 10.42 

19.27 

4 

7.98 

- 2.80 

7.81 

- 3.81 

5 

4.28 

16.29 

5.27 

16.98 

6 

- 0.11 

15.01 

- 1.22 

2.06 



a - 3.20 

o 


1 

- 27.82 

165.92 

- 10.67 

-146.27 

2 

68.31 

- 9.34 

- 50.85 

7.50 

3 

- 6.40 

40.79 

- 18.44 

79.17 

4 

43.78 

- 70.01 

- 24.00 

63.43 

5 

10.39 

28.49 

3.88 

7.05 

6 

- 0.72 

18.02 

- 1.50 

- 0.65 


/Matrix 
I Inversion 
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r int *r ext = 6800 km 


Tnh 1 o 1 


I 
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Internal 

External 



Harm 

Cos 

Sin 

Cos 

Sin 


1 

- 35.09 

22.04 

- 2.92 

- 0.56j 

^Input ' 

2 

22.28 

4.44 

1.86 

- 5.48' 

[coefficients^ 

i 

- 22.95 

26.89 

2.29 

11.16 


4 

18.87 

0.20 

- 2.61 

- 7.75 


5 

7.86 

27.34 

1.74 

6.68 


6 

- 0.70 

17 . 99 

- 1.50 

- 0.61 


- 0.00 






1 

- 35.07 

21.92 

- 1.68 

1.03i 

/^<eeur.slve\ 

2 

22.24 

4.37 

0.24 

0.04| 

l.cast 1 

) 

- 22.88 

26.84 

- 0.06 

0.07' 

lr.«|iiarc-s / 

4 

18.80 

0.20 

0.01 

0.00 

\ / 

5 

7.83 

27.24 

0.00 

0.00 


6 

- 0.68 

17.92 

- 0.00 

0.00 


- 0.46 






1 

- 14.60 

- 0.45 

-424.72 

463. 95| 

^Marrlx \ 

2 

28.26 

12.16 

-551.55 

-629.671 

[inversion/ 

3 

- 34.26 

15.65 

3601.09 

4608.37 

\ f 

4 

22.91 

- 35.98 

6193.85 

64513.81 


5 

9.35 

27.59 

8247.71 

131.75 


6 

- 0.70 

18.02 

- 77.75 

-466.09 


r int “ 6400 

km r ext * 

5 X r Int 






Table 2 








ORIGUmmL page is 






OF POOR QUALITY 





Internal 


External 



Harm 

Cob 

Sin 

Cob 

Sin 


1 

- 35.09 

22.04 

2.92 

0.56 

/input \ 

2 

22.28 

4.44 

1.86 

5.48 

Icoefflclentsl 

J 

- 22.95 

26.89 

2.29 

11.16 

\ f 

4 

18.87 

0.20 

2.61 

7.75 


5 

7.86 

27.34 

1.74 

6.68 


6 

- 0.69 

17.99 

1.50 

0.61 


- 0.22 






1 

- 20.00 

11.64 

1H.‘)/ 

10.75 , 

fKfcurKl v«A 

2 

12.38 

0.40 

12.13 

0.63 

1 I.enHt 1 

\ 

- 10.80 

19.68 

10.45 

18.92 

(Sf|uarfH J 

4 

8.34 

2.79 

7.96 

3.65 

\ f 

5 

4.38 

16.80 

5.19 

16.81 


6 

- 0.07 

13.44 

1.25 

3.76 


1.52 






1 

105.00 

97.15 

148.91 

76.98 

matrix \ 

2 

- 56.05 

54.52 

83.15 

57.60 

VlnverBion# 

3 

32.90 

40.61 

57.87 

4.37 

\ f 

4 

- 12.62 

7.99 

29.46 

12.65 


5 

7.62 

26.21 

1.54 

5.59 


6 

- 0.68 

17.99 

1.52 

0.61 


r int ■ 6400 

r ext ■ 7200 
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Internal 


External 
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Cob 
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Cob 

Sin 

1 

- 35.09 
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- 2.92 

- 0.56 /input 

2 

22.28 

A.AA 

1.86 

- 5.A8 Icouf f Ic lents 

3 

- 22.95 

26.89 

2.29 

11.16 ' 

A 

18.87 

0.20 

- 2.61 

- 7.75 

5 

7.86 

27. 3A 

1.7A 

6.08 

6 

- 0.70 

18.00 

- 1.50 

- 0.61 
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1 

- 32.97 

20.08 

- 9 . 86 

5.87 

/Ui't'urs I ve\ 

2 

21.98 

3.AA 

3.79 

0.5A 1 

[ I.faKt j 

3 

- 22. A5 

27.62 

- 2.17 

2.6A 


A 

18.58 

- 0.19 

1.02 

- 0.12 
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5 

7.90 

27. AO 

0.29 

0.88 


6 

- 0.72 

17.91 

- 0.06 

0.32 



0.18 





1 

- 1A.59 

29.89 

- 67.98 

- 29.73 /Matrix A 

2 

A8.A6 

3A.21 

-151.57 

-167.97 linverslon# 

3 

- 61.68 

36.55 

377.70 

- 6A.65 ' ' 

A 

- 3.96 

- 29. 2A 

A12.28 

500.39 

5 

9.0A 

26.39 

28.25 

- 1A.12 

6 

- 0.68 

18.01 

- 1.95 

1.22 


rint “ 6A00 r ext ■ 2 x rlnt 


Tab 1 «• A 
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Internal 
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Harm 
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Sin 

COH 
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1 

10.05 
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- 5.41 
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2.30 

- 2.91 
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Input 
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0.23 

5 

- 4.90 

- 3.95 

- 0.16 
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6 

11.15 

- 1.27 
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- 0.03 


/uecurKlvf 
I l.faHl 
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A 


V 




/ 


0,00 


8.30 

17.25 

56.03 

6.94 

19.00 

3.46 

83.01 

- 20.25 

42.11 

- 23.39 

-216.52 

129.55 

1.26 

21.81 

57.37 

-332.44 

5.94 

- 4.17 

- 9.39 

0.22 

11.23 

- 1.37 

- 1.51 

3.74 


Matrix 

Inversion 
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Table* •) 
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Sin 

1 

10.05 
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- 5.41 

7.87 /input 


- 3.87 

2.30 

- 2.91 

- A.97 Icoeff Iclents 

3 

20.93 

- 11.23 

- 5.37 

6.39 ' ‘ 

4 

1.21 

3.63 

- 0.58 

- 16.30 

5 

- 5.25 

- 4.04 

8.47 

2.55 

6 

11.23 

- 1.36 

- 1.48 

3.26 


- 1.70 





1 

2.32 

10.41 

2 . 30 

9.13 /ueeurKlve 

2 

5.34 

4.24 

- 4.49 

- 3.691 Least 

3 

13.09 

- 3.41 

0.73 

3. 40 y Squares 

4 

- 2.60 

- 2.25 

- 2.90 

- 19.18 ' 

5 

- 3.11 

- 5.08 

7.87 

3.42 

6 

11.35 

- 1.03 

- 1.33 

3.29 


1.59 





1 

16.42 

22.33 

- 11.06 

2.70 /Matrix 

2 

- 6.66 

5.33 

- 3.90 

- 3.96 linvorslon 

3 

20.02 

- 12.03 

- 5.20 

6.57 ' ' 

4 

1.37 

3.45 

- 0.55 

- 16.33 

5 

- 5.23 

- 4.02 

8.46 

2.55 

6 

11.23 

- 1.36 

- 1.48 

3.26 


r int ■ r ext ■ 6800 
Alt - 6800 + 3400, s In 0 
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